Current resonances in graphene with time dependent potential barriers 
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A method is derived to solve the massless Dirac-Weyl equation describing electron transport in 
a mono-layer of graphene with a scalar potential barrier U(x,i), homogeneous in the {/-direction, 
of arbitrary x- and time dependence. Resonant enhancement of both electron backscattering and 
currents, across and along the barrier, is predicted when the modulation frequencies satisfy certain 
resonance conditions. These conditions resemble those for Shapiro-steps of driven Josephson junc- 
tions. Surprisingly, we find a non-zero {/-component of the current for carriers of zero momentum 
along the j/-axis. 
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Growing interest to graphene, see e.g. Ref. pQ, is stim- 
ulated by many unusual and sometimes counter intu- 
itive properties of this two dimensional material. Indeed, 
graphene supplies charge carriers exhibiting pseudo- 
relativistic dynamics of massless Dirac fermions. As one 
consequence the Klein tunneling phenomenon [2] occurs 
with unit probability through arbitrarily high and thick 
barriers at perpendicular incidence, irrespective of the 
particle energy, in accordance with experiment [3 . The 
question arose how to control the electron motion in 
graphene and hence boosted detailed studies of Dirac 
fermions under the influence of various forms of scalar 
[3HH] or vector [5] potentials. 

Applying a time-dependent laser field to pristine 
graphene opens an alternative and efficient way [TUHT2] to 
control spectrum and transport properties. It was shown 
[lOj that Dirac fermions accross p-n-j unctions can aquire 
an effective mass when driven by a laser field. This re- 
sults in an exponential suppression of chiral tunneling 
even for perpendicular incidence upon the junction, if 
the AC-electric field is directed parallel to the junction, in 
stark contrast to Klein tunneling occuring in the absence 
of the laser field. Actually, time dependent laser fields 
can mimic [12] the influence of any electrostatic graphene 
superlattices on the electron spectrum in graphene. The 
question arises whether and under which conditions time- 
dependent modulations of an electrostatic barrier, where 
energy is not conserved, would affect electron transport 
and generate backscattering. 

In this Letter we answer this question by solving the 
problem for arbitrary space-time dependent scalar po- 
tentials U(x,t). Our solution is based on expanding the 
wave function as a power series with respect to the mo- 
mentum k y parallel to the barrier, and manifests a struc- 
ture of left and right moving waves. All terms appearing 
in the ^-expansion can be calculated analytically, de- 
spite of the fact that each term is described by a partial 
differential equation in (x, t)-space. At k y = (normal 
incidence upon the barrier) we confirm complete Klein 
tunneling for any U(x,t) while for finite k y backscatter- 
ing resonances can occur at certain angles of incidence, 



depending on the modulation frequency of the barrier. 
As a counter intuitive result we find a non-zero and os- 
cillating current j y along the barrier, even at ky — 
for valley polarized fermions. At k y ^ the current 
j v arises also in valley unpolarized situations, it can be 
resonantly amplified and flow in either direction. Inter- 
estingly, j y exhibits a non-zero DC-component at certain 
resonance frequencies, in full analogy to Shapiro-steps of 
driven Josephson junctions. 

At low energies, the honeycomb lattice of graphene en- 
genders two copies, t z = ±1, of Dirac-Weyl Hamiltonians 

m 

H = v F [T z a x p x + 6-ypy] , (1) 

centered about two inequivalent Dirac points ( "valleys" ) 
K and K' at corners of the hexagonal first Brillouin zone 
where electron-hole symmetric bands touch; Pauli ma- 
trices o x ^ ViZ act on two-component spinors representing 
sublattice amplitudes. Carriers near either of the Dirac 
points exhibit opposite Fermion helicities, <j -p/p = ±1 . 
Proposals exist in literature how to valley polarize car- 
riers in graphene, by means of nanoribbons terminated 
by zig-zag edges |14) . by exploiting trigonal warping at 
elevated energies [TH] , or by absorbing magnetic textures 

Smooth electromagnetic or disorder potentials [17j . 
containing negligible fourier components at large wave 
vectors of the order of \K\, will not cause scattering 
between valleys so that calculations can be carried out 
for t z = +1 or t z = —1 separately. Accordingly, time 
dependent potentials U(x,t) should be slowly varying, 
without frequency components that might induce excur- 
sions to energies where the band structure of graphene 
starts deviating from the isotropic cone spectrum, i.e. 
below 0.6 eV Q2]. Including U(x,t), the Dirac equation 
for the wave function tyu (x, y, t) — ^{x, t) exp(ik y y) can 
be written in the form 

{-ir t & U(x,t) ) * + lk v{l ; v l dt [l) 

where from now on we assume v F = 1 and h = 1. 
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This equation has been solved analytically for time- 
independent potentials either by matching [2] of wave 
functions for rectangular barriers, or by the WKB 
method |19j for smooth barriers. Time dependent har- 
monic oscillations have been considered of gate voltages 
on either side of a graphene rectangle [20], of an electric 
field parallel to the barrier |10j or in resonance approxi- 
mation [12], or for some class of time dependent barriers 
U(x,t) at k y = [H]. 

Our goal here is to construct the solution of eq. ^ for 
arbitrary U(x, t) acting at positive times, U(x,t < 0) = 
. From the Ansatz 

* = ( I ) *»,+ + £(**»)" ( \ ) *»,- 

n=0 V ' ' ' »=n V ~" ' ; ; 



n=0 



(3) 

as a power series in k y we derive a recurrence relation for 
the coefficients \P n ,± which obey the inhomogeneous first 
order partial differential equations, 



, . d d 



* n ,±±T,* n _i, T =0, (4) 



with > 0) = 0. Initial conditions can be 

chosen as ^o,±(x,t — 0) = a±(x) = \^A(x,t = 0) ± 
T z ^ B (x,t = 0)]/2, *„> ,± = 0, where *a,*b describe 
electron amplitudes on either of the graphene sublattices. 
The two functions a±(x), providing the initial conditions, 
can be, e.g., a plane wave or a wave packet. Wc underline 
here the general structure of ([3| as a sum of right and 
left ^_ moving waves. Using the standard d'Alembert's 
ratio test, a sufficient criterion for convergence of series 
^ is Ifc^llimn^oo |* n +i,±|/|*n,Tl < 1 for an relevant x 
and t. 

Despite of the fact that Q are partial differential equa- 
tions, we can solve them exactly using the method of 
characteristics [22] . The corresponding result reads 



, ,\ / ,\ —i f dt' U(x^ft±t' ,*') 

(x,t) = a n< ±{x,t)e J a 



(5) 



with ao,± = a±(a; =F t) and 

a„>o,± = Tir z [ dt' m n _ liT (xTt±t',t') 
Jo 

i f* dt" C/(a;=Ft±t",t") 

Together with |3]) the recursive solution for ^ n ,± pro- 
vides the exact wave function ^ to any desired accurancy. 
To zeroth order approximation w.r.t. k y we obtain: 



ip(x, t) = a+(x — t) 



+ a-(x + t) 



1 



-i J* dt' U(x-t+t',t') 



-i J dt' U{x+t—t',t') 



(6) 



The first order corrections w.r.t. k y in (15]) can be written 



Jo 

xg i J o dt" [U(x^t±t",t")-U(x^t±2t'^t",t")] 



so that * 



*+ ( r : J +*_ as in < 3 >> with 

,T, r _L 7 A 1 -i V dt' U{x^t±t' ,t') 



(7) 



(8) 



When k y = and when the wave packet is initially purely 
right moving, a_ = , eq. ([6]) reveals that the electron 
density distribution \a + (x — 1)\ 2 undistortedly continues 
to propagate to the right without reflection: (x, t) = 
for all times t > 0. This proves complete Klein tunneling 
also in the presence of time dependent barriers; wave 
functions ^± acquire only a phase factor by the potential 
at k y = 0. 

As a measurable quantity, we now evaluate the cur- 
rent density in cartesian components, j x = ^*r z a x ^ = 
2#^+ _ 2*1*_ = j Qx + jxx and j v = %*<T y V = 
2ir z (*+*_ - ) = j 0y + ji y . Here, the last equal 

signs refer to zeroth and first order contributions w.r.t. 
k y , respectively, yielding 

jOx = 2(| ao+ | 2 -|a _| 2 ) (9) 



jo y = 4t z |a 0+ a - 1 sin(p + <j> ) (10) 



hx = 4:k y T z nc{a a+ A* 1+ - a -A*_} (11) 



ji y = ik v (\Ai + a* _\ sin(ip - <j>_) - |AJ_a 0+ | sm((p + <£+)) 

(12) 

with ip = Jq[U(x + t - t',t') - U(x - t + t',t')]dt', 
0o = ar g(ao+ a o-)i an(1 4>± = ar g( a o±^izp)- We distin- 
guish two cases: (i) r z -independent contributions jo x and 
ji y which can be observed for valley unpolarized carriers 
and (ii) r z -dependent contributions j± x and jo y where 
detection calls for valley polarization. 

Eqs. ([9]) and ( 10 1 describe the current density at nor- 



mal incidence, k y — 0. Then j x stays unaffected by the 
barrier, irrespective of r 2 which rephrases the above re- 
sult of complete Klein tunneling. Surprisingly, a current 
j y flows perpendicular to the momentum in graphene, 
provided the sample is valley polarized (the total current 
j = j + + j , where j T * originates from states near val- 
ley K (t z = +1) or K' (t z = —1), respectively). This 
current ( 10 1 results from interfering left and right mov- 



ing waves, which both need to have nonzero amplitudes, 
a +a - 7^ 0. 



Eqs. (11) and (12 1 describe corrections to the current 
density at small but finite angles of incidence, k y 7^ 0. 
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Thereby, j\ x exhibits qualitatively similar properties as 
jo y ; in particular it stays nonvanishing at finite valley po- 
larization only. By contrast, the current density ji y now 
exhibits striking current oscillations and current rever- 
sals already in valley unpolarized situations as we show 
in more detail below. 

Next we turn to the question how carriers are re- 
flected by U(x,t). Let's consider an initially right mov- 



ing plane wave, *j 



('.) 



at t = which 



produces a current density jo x — +2 pointing to the 
right. Using equations ([3| and (Jsj) , and assuming small 
k y , the leading contribution to the reflected current den- 
sity ]2x = — fcyl^i-l 2 arises in 0(k 2 ) under the action of 
the barrier at t > and is proportional to \ jax\-, cf. (7|9 1. 
This suggests to employ the ratio 



R(x,t) := -j 2x /j 0x = k 2 y \A^(x,t)\ 2 /2 (13) 

as a measure for the reflectivity at small k y . While the 
quantity R(x,t) evolves in time, together with U(x,t), 
it is independend of r z and, thus, measurable without 
valley polarization. Moreover, we also analyze the time 
averaged reflectivity R(x) :— lim<r->-oo Jq R(x,t)dt/T, 
which can be measured just by means of DC-equipment. 

In the following, two specific examples U(x,t) = 
[/(- 1 ' 2 ) (x,t) are considered. As initial conditions we take 
into account two cases: (i) a superposition of equal 
amplitudes of right and left propagating plane waves, 
a± = cxp(±ik x x) and r 2 = +1 when calculating j 0y ; (ii) 
an incidently right moving wave a + = exp(ik x x) , a_ = 
when calculating R and j\ v for valley unpolarized sys- 
tems. Our first example is U^'(x, t) = Uoxsvuuit . In 
view of (10), we derive for this case 



(14) 



j 0y (x,t)= S m2[k x x + U (±- S -^)] > 
which can be rewritten as a sum 

joy (x, t) = y j J n I I sin I 2k x x H nut 

71 — — GO \ / \ 

(15) 

using Bessel functions J n . This form (151 reveals a pe- 
culiarity at u> = uj n with 



L0 n = \j2U /n , n G N, 



(16) 



similar to Shapiro-steps |23j of a driven Josephson junc- 
tion. As depicted in Fig. la, frequencies u — u> n gen- 
erate periodic oscillations, which, again as in the case 
of Shapiro-steps, induce a nonzero DC-component in the 
current at given x. Modulating the potential with uj =/= uj n 
results in aperiodic oscillations and zero DC-component. 
Similar resonance effects can also be seen in both, the 



reflectivity R (13) and the current j\ y (12). By inserting 
U^(x,t) into eq. (ITJ) we derive 



Ax- = 



(x+t) I ^£ e ~2ik x t' e 2iU (u>t' ' -sin(ut' )) / u) 2 {Yj\ 



1- ' ' ' ' A » A n 




1 V(u ) 



FIG. 1: (Color online) (a) Current jo y (10 1 perpendicular to 
fe versus time for U^(x,i) = Uox smut, k x x = n/8, and 
Uq/uj 2 — 1/tv (blue line) and Uq/oj 2 = 1/2 (red line), as- 
suming a valley-polarized situation t z = 1. For "Shapiro- 
step" conditions (eq. (16l) periodic oscillations can be seen 



(red), while, away from this condition, aperiodic oscillations 
occur (blue), (b) Same as (a) but for potential U^ 2 \x,t) = 
Uo cos(x/L) cosujt with x = ttL/2, k = 0, UoL = 0.1, and 
frequencies ui = (n /2)(1/ L) (red line) and ui — 1/L (blue 
line). Both currents are aperiodic. For the matching con- 
dition uu = 1/L (blue) a considerable enhancement followed 
by a saturation of the amplitude of the current oscillations 
occurs, while away from this res ona nce no enhancement is 
seen versus time, (c) Current ji y ( 12 I for U^(x,t) at x = 0, 
using k — 0. At Shapiro-resonance (Uo/lj 2 = 1/2, blue line) 
pronounced current enhancement occurs, cf. eqs. |l7fT8l ) while 
away from the resonance {Uo/uP = 3/n, red line) no enhance- 
ment is seen, (d) Time-dependent reflectivity R(t), calculated 
by numerical integratio n o f eq. ( |17[ l, blue line, and by us- 
ing the approximation (19 1, red dashed line, near resonance 
for kx/y/Ua = and ur/Uo = 0.49. (e) R{t) at resonance 
(green dashed line, kx/^/Wo = 0,lj 2 /Uo = 1/2) and near 
resonance (blue line, k x /\/Uo = 0,u> 2 /Uo = 0.49). At the 
resonance ^Rit) increases with time ~ i 2 , in agreement with 



equation (19). (f) Time-averaged reflectivity R as a function 
of kx/VUo for u)/\/Uo — 1. Equidistant resonances occ ur a t 
kx/y/Uo = 1 — n/2 (dashed blue vertical lines), cf. eq. (18 1. 



One of t he r esonance peaks is well fitted by the resonance 
equation (201, as shown in dashed green, (g) R as a func- 
tion of the driving frequency uj/^/Uq for kx/VTJo = 0: reso- 
nance peaks are clearly seen at Shapiro-step conditions ( 18 ) 



'2/n, indicated as dashed magenta vertical lines. 
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a ik x (x+t) 



n=—oo 



2Uo\ e i ( 2lJ o/ U) ~' 2k ^~ nu) ) t — l 
Lu 2 J i(2Uo/u> — 2k x — ntu) 



from which we read off a Shapiro-step resonance condi- 
tion 



k x 



Uo 

LU 



n e N 



(18) 



specifying now a directional dependence of the momen- 



tum k. From (17) together with (12) we conclude that 



in valley-unpolarized samples the current j y oc k y par- 
allel to the barrier oscillates as a function of time and 
may take either sign (despite of the fixed k y , see Fig. [Tf;) . 
In addition, the amplitude of these oscillations increases 
with time as the Shapiro-step resonance condition ( [T8| is 
met (compare red and blue curves in Fig. [IJ;) . 

Analogous resonances also show up in both reflectivi- 
ties, R and R. The latter allows experimental observa- 
tion of the here predicted behavior without time-domain 
measurements. Indeed, near the Shapiro-step resonance 



(18) we can keep only one summand in the expansion 



(17), yielding 



2Uo 

2 



sin 2 [([/ 



k x 



j/2)t] 



uj 2 J (Uq/lu — k x — nuj/2) 2 



(19) 



This equation is in a good agreement with numerical in- 
tegration of (17), see Figs. [IJi,e. Averaging (19) with 
respect to time results in 



R — -k 2 J 2 
11 — ^y d n 



2Uo\ 1 

2 I (U /uj -k x - ulj/2) 2 



(20) 



so that the barrier will become intransparent near mo- 
menta k x — k n (see Fig. [ljf), already for small Uq. This 
produces strong anomalies in transport properties at an- 
gles arctan(fc ra /fc y ) of the incidence. Instead of sweeping 
the directions of k one may alternatively sweep u> at fixed 



fc, cf. ( 
in Fig. 
value 



181) ; ensuing resonance peaks are clearly observed 
lg. The constraint R < 1 determins the maximum 



I Ad < 



\Ug - k x - nuj/2\ 
\J n (2U /uj 2 )\ 



(21) 



where second and higher order terms in the expansion 
([3]) can be ignored. 

As a second example, we consider U^{x,t) — 
UqCOs(x/L) coswi to demonstrate how even more in- 
triguing resonance features can arise from the interplay 
between spatial and temporal periodicities. Given again 
the initial condition of left and right moving plane waves 
of equal amplitudes, and assuming valley polarization we 
find 



JOy 



sin I 



x sin 



2kx 
loL - 



4U Q Lsm(f) 



2L 



-t 



LU 2 L 2 - 1 

loL - 



(22) 



sin 



2L 



Now, oscillations of jo y persist even when to 0, since 
the spatial periodicity 2itL of the potential induces a fre- 
quency component v F /L to waves moving at the uniform 
Fermi velocity (restoring here v F ). This reminds of the 
AC- Josephson effect ^2 3J where AC-current oscillations are 
generated by a time-independent voltage. 

On the other hand, if the barrier modulation fre- 
quency uj — > ±u F /L, the argument of the sine in the 



curly brackets (22) varies proportional to t as 2k x x =F 
t Uq sin(x/L) smujt . For small Uq the oscillations of jo y 
thus grow resonantly with time, before they saturate at 
t > 2tt/Uo, cf. Fig. [TJd. We mention the analogy to 
resonant excitations of plasmonic oscillations by spatio- 
temporal mode matching of the incident light with the 
grating period (Wood's anomaly 24J). Similar effects oc- 
cur also for valley unpolarized currents (e.g., j\ y ) and the 
reflectivity R(x,t), but calculations become considerably 
more cumbersome and will be published elsewhere. 

Concluding, we present the analytical solution of the 
Dirac equation for Fcrmions in graphene moving in a 
scalar potential barrier U^Xjt) of arbitrary x- and time- 
dependence. Unit transmission probability, referred to 
as Klein tunneling, is found for normal incidence upon 
the barrier, rendering at most a phase to the wave func- 
tion. On the other hand, under certain angles with re- 
spect to the barrier (k y 7^ 0), we predict strong reflec- 
tion, even for weak potentials. Further, also the cur- 
rent parallel to the barrier, j y , may exhibit oscillations, 
despite of a constant electron momentum k y . The am- 
plitude of these oscillations grows linearly in time when 
U(x,t) meets certain resonance frequencies. In valley- 
polarized samples j y does not vanish even for zero mo- 
mentum parallel to the barrier {k y = 0), provided left 
and right moving waves both interfere with finite ampli- 
tudes. For graphene nanostructures driven by oscillating 
potentials, the predicted resonances in current and re- 
flectivity can be seen, for example, in electron transport 
properties (e.g., in AC and DC electrical conductivity) 
through suitably arranged quantum point contacts. The 
new non-stationary phenomena in graphene calculated 
here within the single-particle approximation can pro- 
mote development of a more elaborated many-electron 
non-stationary theory of AC-driven graphene nanostruc- 
tures which is crucial for future graphene-based electron- 
ics. 
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